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INATIONS OF RANDOM LOADS 


“Many physical structures are threatened by combinations 
of loads of varying ragnitudes from various sources. For 
anstance, bridges, piers, dams, and buildings can experience 
loads frem wind, snow, ice, tides, earthauakes, and so 
forth. In many instances the total load or stress exferi- 
encec Fy a structure varies in time in an apparently randon 
fashicn. Certain components of the loads vary rather 
Slowly; others occur more nearly as impulses, such as those 
associated with winds or earthquakes. The problem is to 
design a structure to withstand the superposition of random 
loais from many sorrces with at least an approximately 
understocd probability. 

The purpose of tkis thesis is to describe and investi- 
gate certain simple but somewhat realistic probakilistic 
loal models for use in design, and perhars safety assessment 
of structures. In this thesis we confine attention tc the 
superposition of just two load types; shock loads, and 


constant loads. For example, winds gusts, flash floois, and 
earthquakes have varying magnitudes and have relatively 
short durations in cecmparison to the times between their 
eccurrences. These will be modelled as instantaneous shock 
loads. Cn the other karnd, snow, ice, or water accumulation, 
or even the presence of slowly moving vehicles present loads 
that remain nearly ccnstant in time, occasionally changing 
to new levels. These will be modeled as constant loais that 
change inifreguently. Throughout this investigation it will 
te assumed that the effective stress exertel by several 
types of loads acting Simultaneously can le expressed as a 
linear combination cf the component loais, the load 


conporents being treated as stochastic processes. 


2 


The times between changes in magnitude of the constant 


load process are independently and identically distrikuteda 
with distribution function H. The successive magnitudes of 
the constant load are independently and identically distrib- 
uted with distribkuticn function F. Given the constant load 
process, the shock load process is a compourd Poisson 
process. The successive shock load magnitudes are iniepen- 
dently and identically distributed with distribution func- 


tion G. The conditicnal rate of arrival of shocks, given 


.that tle magnitude cf the constant load process is x, is 
MX) - 


Let X(t) be the magnitude of constant load precess at 
time t; anc Y(t) be the magnitude of shock load process at 
time t; Z(t) =X(t)+Y(t) is the superposition of the two loads 
at time t, ard H(t) =SUE Z(s), the maximum load combination 
in (C,tj, see figure 1.1. 

Let T,y={inf t20:; 2Z(t)>x } die, the first-passage time 
for the loac combination to exceed a stress level x. aa 
will represent the time to faiiure of a structure whose 
strergth is x, and is subjected to a stress history 
{Z(t); t20}. 

Gaver ard Jacobs (1981) studied the above inodel ror the 
case in which the distribution H is exponential and the 
shock load rate A’ is a constant independent of the shock 
ioad process. Related load combination models that have 
Leen studied include Peir and Cornell (1973) i Ref. 1) Wen 
(1977) [Ref. 2] Pearce and Yen (13983) [Ref. 3]. 

Asytftotic results which appear in Gaver and Jacobs 
(1924) and Jacobs (1984) will te described in chapter 2 for 
the distrikution of Tx aS xo ana for the tail of the 
cistribution of T, for finite x. 


In chapter 3, a Simulation program for the load 


conbination model will be described. 
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The Load Combination Process. 


Simulation study of the use cf 


of 


chapter 2 to approximate 


of Ty will be described. 
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II. ASYMPTOTIC RESULTS FOR THE DISTRIBUTION OF 


A. THE DISTRIBUTION CF Ty FOR LARGE X 


This section summarizes some of the results founé 


Gaver anc Jacobs (1981) fRef. 4] for the model in which 


in 


the 


Eistrikbution H is exponential(a) and the shock load arrival 


rate A is constant independent of the constant load frocess. 


First, we consider the distribution function of the 
maximum load combination to occur juring £9,ti, M(t). 
B(t) =P {# (t) $x} 
x 
=PE"(t) Sx, Tr >t} +P (M(t) <x,T) St} eee eee (2-1) 


where T,; is the arrival time of the first shock. 
A renewal argument yields; 

(M(t) <x,T) >ty= e*fTexp(-uet or-4) pcdy) 

2M (t) Sx_% <t} = [Fredy] * exp(-av Gerd) Hxct-v) Pedd) 
where G(x) =1-G(x). 


Next take Laplace transforms with respect to t of (2-1). 


Hig) = ff eft Hee dt 


de eee 0k 2 ate ins eee te (2-2) 


where 


MU =f (E4a 4 MG Oey) TF (Ay) cree ne ene ene nen nnn nee (2-3) 


Fence 


ae ( 
BAB Y= TO ay) 


Tt seems to be difficult te invert Hy (E) for any inter- 


estince choice of the distributiors F(x) and G(x). However, 


uselul information can still Le gleaned from equation (2-4). 


First, from the jlefirition of Tx, 


EQM(t) Sx}=P(Tr>t) 0 re eee nee ----- (2-5) 
Thus, Hler= fe St p(t >t) dt  -------------------- (2-6) 


Tet § 70, in eguation (2-6) to find that 


11 


w(x) =E(%) =H (0)= AAS) - = Bae(2-7) 
The next result is the limiting property for the first 
time the load combination process exceeds a given level x. 
Tne limiting distribution of Ty is exponential in the sense 
that 
Tiere (sy Bese 0 HERE (2-8) 
where x, =inf (t;F*G(t)=1} 
The precf of this result for a more general model that 
includes that of chapter 2 can be found in Jacobs (1984) 


[Ref. 5}. 


B. THE TAIL OF THE DISTRIBUTION OF Ty FOR FINITE X 


In this section we will summarize results for the asymp- 
totic behavior of P(I,>t) for large t and fixed x as tow 
for the model in which tke distribution H is exponential and 
the shock load arrival rate is constant. More details and 
results ‘fcr a more general model that includes that of 
chapter 2 can Le fund in Jacobs (1984). 

From the equations (2-1) and (2-5) 

EB (Ty >t) = @-rt{* Ecdy) MEH 


+ rf Pers ds [re cdy ev MEOWS Petr ts) wooo (2-9) 
iet 
L(4s)=f" ody) rneTATHROD)S ds wee (2-10) 
I (t) = i Fcdy) OTIS t- e tMaury) =) SSeS S= (2-11) 

r» 
L(o)= {* Fcdy) NFaeorgy ------- (2-12) 
The renewal equation (2-9) gives; 

P(T, >t) =g(t) + [Freds) Hy (t-s) wae -H~ 12-13; 
where g(t)=[tpcdy) eTrAtmacxp)t ir oat (2-14) 


Following the development in Feller (1971; page 376) 
{Ref. 61 it will be assumed that there exists a number K 


such tkat 
te JP eX? Ledy) nee en nn anne n------- (2-15) 


eo AY ea ere as 


pw ¥ --s — w Sa ee eee ek oo ee a a ead Real, Sa ia in Renn ata iat in Ee eee ares ee 


This roct is unique ani since the distribution of I is 


T . defective, is positive. 

Let 
L¥(dsy= enS Leds) eee een eee ee (2-16) 
- g*(ty= ext git) we nnn eee -------- (2-17) 
a PRetyoty= erPtH, (t) teen nee (2~ 18) 
: Then 

P* (Ty >ty=g%(t)+ fte*(t,>(t-s)}L* (ds) --------- (2-19) 


It is assumed that the assumption for the key renewal 


theorem tolds. Aprlication of the key renewal theoren 


yields 
* KOOt m athe COPE AS: Neuer ee cunteSe a 
Lime P(T, >t) ra [9 (t)at (2-20) 


G4 (tydt= ['F(ay) ("eX erradow)t at 


=f F (dy) Ata Gle-§) - KX) 


Rewriting the equation (2-15), the equation determining 
Kis ' 
x - 
r= [2 eX5 fF ecayy ae LAtMEUWIS de 
A 


Pe eget ES ( 2= 22) 
and 
pts (Pel reda) exert neler Uax- jt dt 
= JX Edy) TAS “Senarz Kop (2-23) 
the key renewal theorem implies; 
Lin e*t 9(7, >t) = ae Ho meee oe (2-24) 


ae No Fay) ———_*_______ 
Examples: [A+ M @(x-4) - KO) 

It appears to he difficult to solve for K analytically 
in general. We will discuss the numerical computation of XK 
for two examples. In both examples, A=AM=1 and the 
Gistribution of the shock and constart load magnitudes F and 
G are of the form; 

F (x) = € 7% 

Gixyeces 


5 CT a eS er OT ee ee ee 


In example 1, a=b=1 and in example 2, a/b=2. 


To compute the kK for F(x)= e°* ,G (x)= «™ . 
let 
PAN 
me. I Se ee EE SE A CNT = 
£G) (Fh ae exw (2-25) 


where A= =1. 

Fer model 1, 

f£ (x)= Trae AU 20- 079-0 Ln 1 4078) = Irn e*)}} -- (2-26) 

for K<1+e* and x ¥1. 

For K=1, E(x) has a singular point, so by using the 
L’HOFITAL RULS [ Ref. 7] 

sin f(x)= Lin g'{x) 7h {k) 
where g9(&) is the numerator and h(x) is the denominator 
cf the function F(x). 
For x =1, 
£(1)=% (&*-€*) 2205:(9-97) 
For model 2, 
£ (x)= TARE {uinkte Dy" (ume re Hy” ze ((1- Kt) 


+e als (1-K+ ¢7O) — Bn ie" + e™)}} sso> (2-28) 


for x <1+e™ and K #1. 
ror K=1 £(K) has alse singular point, so by using the same 
rethcd as used in model 1; 
fork =1 
£(1)= gap (e*-e™) ---- (2-29) 
Lifferentiation of the ecuation (2-25) gives 


rg = {[* 9-94 et Se he sce Sees e ae - 


Eyuation (2-30) is greater than G for Osx<1+te@andx*1. So 
the function f(x) is a monotone increasing function in the 
interval 0sx<1+e™% - Using this result, we compute the k 
from the equation (2-22) by applying the bisection method 
[Ref. 8]. 

Using the computed k we compute the constant ¥* such that; 


nae tes ©: eee 
: J. Pete C A+ GU-4) -KJZ 


is the censtant of equation (2-24). 


14 


ewe te 


rer case a=b=1, 


. rie Tams {lt 40%) Zen -K 40-4) - ee ] = 
~[Le-%(2-K)-2e7* In (eE™(2-K)) — mas) -- (2-31) 


For case a/b=2, 


UeKseby® (u-wetseX® 5. 
{Li-weteby® (ume tte 5 ob (-xy(1- 0-6) 


4. 
r (1-K)* 


43.0720 tn( 1 e407) In (UK ee OME Oe aes g pal 2> 32) 


ITI. COMPUTER SIMULATION 


this section describes the computer simulation model. 


we denote the constant load arrival rate bya”, the shock 


load arrival rate byy, the distribution of constant load 


Magnitude by F(x), and the distribution of shock load magni- 


tude by G(x). The simulation model will be used to study 


the quality of the approximations of the distribution of 


by the asymptotic results described in chapter 2. 


The following describes the simulation algorithko. The 


StreS€ i€VEIS are X, <X,<Kz oe ee eee KX ye 

1, Set T=0, N,=0. 

2, Set V=0, generate a constant load magnitude,C, and 
duration interval Tc. Find the largest index i>N, 
such that x; <C. 

Put N, equal to that index and 
set 

Ty, =T for N,<is<N,, 
set No=N,. 

3, Generate a shock load magnitude,S, and a time until 
the shock occurs,Ts. 

4, If the shock interval does not exceed the excess life 
of the constant load duration,Tc, fine the largest 
index i>N, such that x; $C+S. 

Let N, be that index. 
Set Ty, =T+V+Ts for N, <i<N,. 
Set N, =N, ,V=V+Ts, Tc=Tc-Ts. Go to step 3. 

5, rf the shock interval exceeds the excess life time 
of the constant load duration, 


set T=T+Tc, go to step 2. 
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eGR gae ROL ete tae ——— 
> 

' 

. 

\ 

. 

e 

. One realization cf the simulation is completed when one 

{ : Ty is computed for each level x, ,1Sisn. Zack sirulation = 


consisted of 5000 replications. Further details are 


eee 


(a Figure 3.1 Generating T. EN 

| 

: included in chapter 5. 

- The key subprograms are following. Other programs will al 

af be explained by reading from the output results. as 

FAXDC ; generate the Tx. ree 

7 CC’PUT 3; compute the P(Tx=0), true P(Tx=0), sample nean, Poe 

true mean, standard deviation, coefficient of ae 

° variance. 

| Sample probability 

: Ei We-se ne is obtained from } 
the sample. z 

oe True protability 

T (Tx =0) =1-P(Tx>0) =1-2 (M(0) Sx) =1- fF (ay) =F(t) is ne 
oktained from the equation (2-5). ie 


The sample mean 
= n 
le x aa Tx: and 
feat 


we oe pie ak ar ere ee SE Seg etre Sree Sad han eh Se os ee eit ST 8 lt aa sages We Sota ee Ses ae, Sra ee ag er eee ee ae ed 

* * _ . . - . . ~ . - ° . 7 . . ee eS . . . . . . . . ~ . + - . . . - . . . . . . 
aed a . a se MS eg a ie ee a> oe. - aw “a etn fe a . a, 1S, 2 te Me Te "Ce Te ae aes = moe eee ~ . Pears . Per tat . i ae ‘ « 
Pe ML OL a a Re Rae ee ee a ee Mam at et de ee ee ete Me eS tee ee he Bae te Be te 


the true mean E(T,) 1S obtained in special 


= 


cases from the equation (2-7). 
See details in Gaver and Jacobs (1981). 
Sample standard deviation 
i = Hor & (te-*) is obtained from the sample, 
and 
Sample coefficient of variance 
Ctv) = —-- is obtained from the samfle 
mean and sample standard deviation. 
+! CCMP 3; compute tke sample guantiles, and the quantiles 
: cf an exponential distribution with mean E(x). 
A sample quantile is 
%, =(kxsanple size) th order statistic of sarple anc 
e fd = ~E (te) Im( 1-4) is the quantile of exponential 
distribution with mean E(Tx). 
ccm compute the quantile of the sample conditional 
ution with mean 2(Ty)/P(Ty>0). 
The sample guantile of the conditional distribut- 
ion 
Gc=L Ax (number of positive Ty) ]th order statistic 
is obtained for sample cf positive Ty. 
The conditional quantile of the exponential dist- 
ribution with mean E(TIx)/P(Tx>0) is 


Ett) 


hue= ~ Peayaqy TRU) - 


SK 3; compute the K such caked 
K=( K(x) 2 =f" ecdy) <— 
bisection method. 


xt gia ParTET ee, by applying the 


INT ; pes the constant y* such that 


x yes Reed _— PN 3 -——— 
‘e : #) SOT Ser Ee 
by egnevieaeine the K found in SK. 


18 


ai stribution given Tx>0 and the exponential distrib- 


Pg tia 


The following cases were simulated. 


g 

q Listribution cas€é Av A) F (x) G (x) 

. a,exronential 1, 1 1 1-e™ 1-e% 
a (a=b=1) Z, 1 2 " " 

a; 2 1 " " 
b,exponential 1, 1 1 1-27" 1-e%* 
| (a=2, b=1) - 1 2 " " 

¢ ay 2 1 " " 


| c,expenential 1, 1 1/e 1-e7* 1-e™% 

: with varying 25 1 es " " 

| arrival rate 3, 1 et " " 

° 

d,Fareto case 3 1 1-3; ~—Cti«CAT-e™® 
2 1 eras 1-e™% 
3, 1 i a Sabre 5a 


is These different settings will be referred to as case (A-1), 
(A-2), (A-3), (8-1), (Br2), (3-3), (C-1), (C-2), (C-3), 
(0-1), (2-2), (D-3) respectively throughout this paper. for 
each akove case, we ccmpare the asymptotic result (2-8) and 
the experimental results. 

The next setting is to examine the exponential approxi- 
mation of the tail cf the distribution of T, for finite x. 
The cases of identical exponential distributions, 

F(x) =G(x)= e” and different exponential Cistributions, 


F (x) =e mde 


G (x) =e where a/b=2 

were considered. 

The following conditicns are considered. 
e, dr =A = F (x) =5 (x) =e™ 


rae NBM =I P(x)=e"2% = ,G(x)=e™* 


iS a eS ea ek at RES se Mla ae ah ey eRe Ses er gage ee, 
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These different settings will ke referred to as case (E) 


case (F) 
constant 


results 


respectively. For each case, compute the K, 


and 


and . - 


v* and compare the Juantiles of the asymptotic 


(2-24) to the simulated data for each level of x. 
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As indicated in chapter 2 the exact distribution of 7, 


is difficult to obtain in gereral. In this chapter simula- 
tion will he used to study the accuracy of the approxima- 
tions of the distribution of TY, by the exponential 


distributions suggested by the asymptotic results of chapter 


A. THE DISTRIBUTION CF Tx FOR INCREASING X 


The limiting result (2-8) indicated that the distribu- 
tion of Tx approaches the exponential as level x increases. 
This result suggests that the Gistribution of Tx can be 
approxinated by an exponential distribution at least for 
large x. In this section this exponential approximation will 


ke studied via a simulation. 


1. Constant and Shock Load Magnitudes Have "dentical 


In the first collection of three motels to be 
considered F(x) =€(x) =e. The times between constant load 
changes are exponential with parameterAand the sheck ioad 
arrival rate « is a constant independent of the constant 
loai process. 

Tor these cases it is possible to derive an analyt- 
ical expression for E({Tx). 

F tg )= $e [rereen- ater sete hte) | 
So at ed « 


In the Tables, P(Ty =)) refers to the simulated 


sample probability that Ty =0; xX-BAR gives the simulated 


sanyle mear of Tx, 3; ST-DEV is the simulated standard 
Geviation of T,; VAR(X-3B) is the variance of the sample 
mean; CCEF-V is the sample standard deviation divided by the 
sample wean. T P(7,=0) is the theoretical prohability that 
Tx=0, (x); T x-bar is the theoretical mean of Tx compute? 
irom the equation (4-1). 

The simulated coefficient of variation in all three 
caSes decreases as the level x increases becoming quite 
close to the theoretical exponential value of 1 when x=5. 

To compare the simulated distribution of Tx to the 

pproximatiny exponential distribution, quantiles fror the 
Simulated fata were computed and compared to the corre- 
sponecing yuantiles cf an exponential distribution having 
theoretical mean (4-1). The Simulatei yuantiles were 
computed as described in cha,ter 3. The exponential 
oa-juantile is given ty 
Gy =7a-lr (1-4) 

were a is the mean cf the exvonertial. 

For each level of x, the first row in the tables 
(hei), (A-2), (A.3) for the guantiles of the distributicn of 
Tx gives the quantiles of the simulated jiata and tke second 
Tow gives the corresponding exponential guantiles for an 
exponential distribution having theoretical mean (4-1). 

Cne wav the Gistribution of Tx differs from an exfo- 


nential is that it has an atom at 0. In particular, 


x 


2 (Ty =0) =F (x) =e" 
for the case considered here. 

Cuantiles were usei to compare the simulated condi- 
tional distribution of Ty given T, >0 to an exponential 
Gistrifution having mean F(T. )/P(Tx>0). The Simulatei guan- 
tiles fer the conditional distribution were comrute? as 
Gescribed in chapter 2. These gantiles were compared to 
those of an exponential distribution having theoretical nean 
E(Tx) /P (Tr >C). 
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In the tables (A.1), 


tiles of the conditicnal distribution for each level x, the 


(A.2), and (A.3) for the quan- 


first row shows the simulated quantile and the second row 
the corresponding approximating exponential guantile. 

The quantiles of the simulatei conditional distribu- 
tion are much closer to their exponential approximation than 
the quantiles for the unconditional distribution to their 
exponential approximatior. Fowever the gquartiies for the 
unconditional distrifution get closer to those of their 
approximating exponential as x increases. 

Comparing tables (A.1), (A.2), and (A.3), it appears 
that increasing the arrival rate of shock or the rate of 
change of the constant load magnitude decreases the quan- 
tiles. The change of arrival rate of shocks appears to have 


the greater effect. 


2. _lonstant and Shock Load Magnitudges Have Different 
Exponential Listributions. 
In the next three cases, the distribution of the 


shock load magnitudes is exponential with parameter 1 and 
the Gistrikution of the constant load magnitudes is exponen- 
tial with mean 9.5. The other model assumptions in tatles 
Re1%, (B.2), (R.3) correspond to those in tables (A.1), 
(A-2), and (A.3). 

As before, for each level x, the first row ir the 
guantile table gives the simulated guantile and the second 
row gives the approximating exponential guantile using the 
exponential distribution with theoretical mean. Comparing 
the guantiles of in table (B.1) with (A. 1) (respectively 

B28 with (A.2) and (2.3) with (A.3)), it is seen that 
Gecreasing the mean constant load magnitudes has increased 
the cuantiles. Turther, it appears that the convergence of 
the distribution of Tx to an exponential is faster in the 


case of smailer mean constant load magnitude. 


23 


. 


The exponential approximation to the conditional 


Gistribution of Ty given T,>0 has theoretical mean 

E(T, )/P (Tx >0). Once again for small levels <x,it appears 
that the exponential approximation to the conditional 
Cistributior of is better than that for the unconditional 


one. 


3. _Arrival Rat 


In the next 3 cases, constant load magnitudes and 
stock load magnitudes are exponential with mean 1. Constant 
aoad magnitudes change according to a Poisson process with 
rate 1. Given the constant load magnitude at time t is C, 
the prokakility a shock load will arrive in the time 
interval [t,tth] is m«(C) hto(h). 


In case (C-1) the ccnditional shock load arrival 


rate is m(c)'=Cc. Comparison with table (A.1) indicates that 
cond@itional shock arrivai rate wc)! =C has the effect of 
cecreasing the quantiles of fT x. Further the exponential 


approximation to the quantiles of Tx, is not as good as in 
case (A-1). 
The approximating exponential distribution to the 
conditional distribution of Tx given Ty>0 has a mean of 
T(t )/P (Tx >0). The quantiles 9f the exponential approxima- 
tion to the guantiles of the conditional distributicn afpear 
to be closer than those for the unconditional distrikution. 
In case (C-Z) snock ivads arrive with conditional 
arrival rate «{c)' =exp(c). Corparing the table (C.1) ani 
(C.2), it is seen that the quantiles of the case w{C)’ =C are 
less thar those for the case at(Cy=exp (Cc). This is Fecause 
interarrival times tends to be larjer for the case 
AM(C)?T =exp(C). Compariny tables (A.1) and (C.2) indicates 
that the exponential approximation to the quantiles cf “, is 


not as good as in the case M=1. 
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In case (C-3) the conditional shock arrival rate is 


M(CY' =exp (-C). Comparing tables (C.2) and (C.3 indicates 


that the quantiles fcr the case 

A(C)’ =exp(-C) are less than those for the case y(C)' =exp (C) 
ard furthermore the exponential approximation of Tx aprears 
to Fe clese for the case (C-3) than that of for the case 
(C-?). Ir all of the cases the simulated mean was used as 


the pa er in e exponenti ximations. 
the parameter th I tial approximat s 


In the next three cases, the times between constant 
loai changes of magnitude are exponential with mean 1. Shock 
loads arrive according to a Poisson process with rate 2. In 
case (9-1) the constant load magnitudes have a Pareto 
distribution with narameter =1 and the shock load magni- 
tudes are exponentiaily distributed with mear 1. Comparison 
with table (A.2) indicates that the Pareto constant loa 
magnitude has the effect of decreasing the yuantiles of Tx. 
Further the exponential approximation to the quantiles of i 
is not jood in the Pareto case; (the approximating expenen- 
tial has a mean of tte simulated sample mean). 

"he approximating exponential cCistribution of %& 
given | >0 to the conditional distribution has a mean of 
F(T, )/P (Tx >9). Once again the gquantiies of exponential 
a, proximation to the guantiles or the conditional cistriLu- 
tion appear to he cleser than those for the unconditional 
distribution. 

ar case (D-2) the constant load magnitudes have a 
Tareto Jistribution with parameter 4=2. Comparing takies 

Del) anc (9.2) at is seen that the jyuantiles of the sinu- 
late distribution cf Tx for the case &=2 are larger than 
those for the case d=1. Comparing tables (N.2) and (A.2) 
indicate that the exzyonential approximation to the guantiles 


of Tx is not as geod for the Pareto case as for the 


exponential case. 


In case (D-32) the constant load magnitudes have 
Pareto distribution with parameter «=1 and the shock load 
Maghitudes have Pareto distribution with a~=1. In all three 
cases the Simulated mean was used in the approximating guan- 
tiles. 


B. TRE EXPONENTIAL TAIL OF THE DISTRIBUTION OF Ty FOR 
FINITE X 


In this section simulaticn will te used to study the 
exporential approximation 


suggested by the asyrptotic result 


Tin Péme) =, 
t>0 ext r 


This is an apvroximaticn for P(T, >t) for fixed finite x; it 
Should become more accurate as t becomes large. 

Two caseS are considered. In both cases shock loads 
arrive accorjiing toa Poisson process with rate 1 and 
constant load magnitudes change at tne times of arrival of a 
Poisscr frocess with rate 1. The shock load magnitudes have 
exponential distribution with mean 1. In case (E£), the 
constant loai magnitude distribution is exponential with 
mean 1; In case (F) it is exponential with mean 0.5, 

Lescription of the computation of « and for these two 
cases can Fe found in chapters 2,3. It can be seen from the 
compute] value of K(x) and v(x) appearing in tables, E.1) 
and (F.1) that X(x) is approaching the value 1/E(Tx) (where 
Z(t) is the theoreticai mean) as x becomes larger and for 
tue cases computed K(X) $ Fe - Further , ¥(x) is approaching 
1 as x becomes larger. 

To study the exponential anproximation to the distribu- 
tion of T , quantiles fron the simiiated jata were conputel. 
These can be fourd in tables (2.2) ard (7.2). for each level 
x, the first row presents the simulated juantile; the second 
row cives the approximating e-dauantile computed by 
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The third row shows the approximating X-qguartile computed by 
A 2 1 Na (i ) 
where the theoretical mean is used for E(T,)- 

Comparing the guantiles in tatles (£.2) and (F.2) it can 
te seen that the approximating quantiles (ZL) are close to 
the Simulated ones even for &x.40. The approximating quan- 
tiles computed in (yg) do less well but improve as x becomes 
larger. The two approximating quantiles approach one ancther 
as x becomes larger, as expected. Approximating guantile (z) 
has the advantage,however, of beiny easier to compute. 

Table (G.1) presents 90 * confidence intervals for 
dA7~guantiles with 420.9 for case (E). These confidence inter- 
vals were computed using a large sample approximation, see 
Conover (1980; page 111-114) [Ref. 9]. 

Comparisen with table (E.2) suggests that the approximating 
guantiles approximates the simulated ones well for these 


values. 


C. CCNCLUSIONS 


Approximating guantile L is aiways better than approxi- 


mating guantiley. Approximation L approximates q, well for 


(aS SMail as O.49. However ~ is more difficult to compute 


thank . 
The performance of approximating quantile Z can be 
improvec by changing it to 


= _ _ EC) 


Fe PCT«20) ee 


and using it to approximate the quantiles of the conditional 


Gistribution of "y given Ty >0; for the models considered 


P(™ =0) is easy to cempute being equal to F(x). 
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TABLE A.2 
Identical Exponential( case A-2 ) 
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Different Exponential( case B-1 ) 
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Different Exponential( case B-2 ) 
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V. SIMULATION ALGORITHM 


A. VARIABLE DEFINITICN 


Tx; ={inf t20 ; 2(t)2xi } 

where Z(t)=superpcsition of constant load and shock loai 
Nrow=maximum number of level X; 5 X,<X,<%,<- eee ee ee ee eK 
Ncol=maximum number of repetition of generating T 
Lce=constant load magnitude 

Ls=sheck load magnitude 

Tc=constant load interval 

Ts=sheck load interval 

Typc=ccunting the number of repetition 


L=counting the maximum level of X; 


E. AIGORITHM 


Inputs;parameterS (rA,vuae a, by, seed numbers ) 
Initialize the stress level Xj 3 X, <X,<Xy<-0 eee eee MK 
ipc=0 

Fepeat 


set T.=0 


Repeat 

generate ccnstant load;Le 
generate constant interval;Tc 
7Tind L such that L={sup i;Le2x; } 
Set Tx; = 70 for 1S$is<L 

M=L 

Ci=Te 

A, generate the shock loaé;ls 

generate the shock interval;Ts 


, If Ts>Ci, then 


ty 


Sey 


eee a le PS OE oe we Pee Pome etc 8, heh wl ese Dm oa aw Toe Seas SU SR RN Se ae oe ES : . * 
- et ea a a a nt ete ee ee ee ae ba te tee tt te Pe te mmf a* 2° a 


TO=TO+Ic 

Ts=Ts-Ci 

TTs=0 

generate constant loadjLlc 


generate constant interval;Tc 


Find L such that L={sup i; Lc>Xi} 
> 


If L>M, then 
Set Tx; =T for M<isL 
M=L 
Else ccntinue 
Ci=Te 
Go to B 
Find I such that L={sup isks+Le2n } 


TTs=Tis+Ts 
Set Tx; =Te+TTs for M<i<L 
Ci=Ci-Ts 
Go to A 
Jntil (L=Nrow) 
Ipce=Ipc+1 
Jntil (Ipc=Ncol) 


End Aigorithn 
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COMPUTER PROGRAM 
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